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We establish a general framework to explore parametric statistics of individual energy levels 
in disordered and chaotic quantum systems of unitary symmetry. The method is applied to the 
calculation of the universal intra-level parametric velocity correlation function and the distribution 
of level shifts under the influence of an arbitrary external perturbation. 
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Ensembles of random Hamiltonians are used frequently 
to model properties of diverse physical systems. Al- 
though eigenfunction statistics can play an important 
role in some applications, typically one is interested in 
the statistics of the spectra {ei} of such Hamiltonians. 
Among all the possible statistical ensembles |l| a spe- 
cial role is played by the three invariant Dyson distri- 
butions 2J. Characterized by the assumption that the 
distribution is invariant under, respectively, orthogonal, 
unitary, and symplectic rotations of the basis, the ran- 
dom matrix theory (RMT) has proved to be remarkably 
successful in modeling physical systems ranging from nu- 
clear spectra and mesoscopic quantum dots 3 to in- 
dividual chaotic quantum structures 

An important class of problems arise when individual 
members H of the statistical ensemble undergo paramet- 
ric evolution according to the rule H — > H' = H + XV, 
where V is a fixed matrix, and X is the strength of the 
perturbation. Instead of the random variables ej, one is 
now confronted with random functions £i(X) of the ex- 
ternal parameter X. Once cast in terms of the rescaled 
variable x = X^J ((dx^i) 2 ) (all energies being measured 
in the units of the mean level spacing A), it has been 
argued .6| that, for a generic perturbation V (see below), 
the statistical properties of the entire random functions 
6i(x) exhibit the same degree of universality as that of the 
parent Hamiltonian H. As well as mesoscopic and chaotic 
quantum structures, universality of the random functions 
&i(x) finds application to a variety of physical systems 
including step configurations on vicinal surfaces pj , non- 
intersecting random walkers in one dimension Q , and the 
world lines of one-dimensional fermions [|| . 

Beginning with the seminal work of Dyson Q, the 
statistical properties of the random functions jejfa)} 
have been the subject of numerous investigations [l(j 
(for a review, see, e.g., Ref. 0). To date, exact ana- 
lytic expressions have been obtained for the distribution 
of "local" (in the parameter space) properties such as 
"level velocities" d x ei(x) and "level curvatures" 

d 2 €i(x) At the same time, parametric correlations 

between the sets {ei(x)} and {ei(x + x)} have been ex- 
plored using field theoretic techniques [TJ]. As well 
as establishing the range of universality, explicit expres- 
sions for the parametric correlator of the density of states 
(DoS) and the related level velocity correlation function 
c(u,x) = (Y,ij9x(-i{x)d s ej(x + x)S(e - ei(x))S(e + u - 



ej(x + x))) have been inferred. 

When RMT is applied to many-fermion systems, an 
important distinction arises between two classes of corre- 
lation functions: Employing the terminology somewhat 
loosely, the function c(u>, x) can be termed "grand canon- 
ical" in the sense that the level ej (x + x) need not be a 
parametric "descendant" of the level ei(x). Indeed, the 
definition of c(u>, x) tracks correlations between €i(x) and 
a "descendant" of any other level. However, often one is 
interested in the parametric evolution of the Fermi level 
or a low-lying excitation in systems with a fixed number 
of fermions. Inasmuch as such a level can be interpreted 
as a single particle level of some effective random many- 
body Hamiltonian, the relevant objects are the "canoni- 
cal" correlation functions, exemplified by the intra-level 
velocity correlation function 

c(x) = (d s e t {x)d s ei(x + x)) . (1) 

A different perspective is provided by the distribution of 
single-level shifts 

p(uj, x) = (S (ei (x + x) - ei (x) - u>)) . (2) 

Apart from their intrinsic interest, it has been argued 
that the functions c(x) and p(u>, x) describe parametric 
correlations of resonant conductance peaks of quantum 
dots driven into the Coulomb Blockade regime and per- 
turbed by a magnetic field or an external gate potential. 
Similarly, once x is identified with magnetic flux, c(x) co- 
incides with the ensemble-averaged correlation function 
of single-level persistent currents • 

To date, studies of parametric correlations of individ- 
ual energy levels have been limited to numerical investi- 
gations. Despite its affinity with c(w,x), the intra-level 
correlation function c{x) (and its canonical counterparts) 
belong to a different class of objects. Its analysis presents 
technical difficulties which are, in part, similar to the 
challenges encountered in the calculation of the level 
spacing distribution in non-parametric random matrix 
ensembles. The latter are known to engage DoS correla- 
tions which go beyond the two-point averages presently 
accessible by field theoretic techniques, and lead to ex- 
pressions in terms of Fredholm determinants with inte- 
grable kernels 0. 

The aim of the present letter is to formulate a general 
framework to explore parametric correlations of individ- 
ual energy levels. In particular, for the unitary random 
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matrix ensemble, we will show that 
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where z±(</>) = 1 + e ±4 * and, adopting the shorthand 
notation 1 to represent the Dirac (5-function, 

Z(cj,x;4>) = det (t-K.(<f>f) . (5) 
The operator kernel K, ((f>) has matrix elements 
^(A,,;^) = ^^(A)^) f(A; ^ F(/i ^ ) , (6a) 

47TI A — jti 

where 

f (A; 0) = fd^^M- (z + z-)V-\\) 

(6b) 

Here the integral is understood in the sense of the Cauchy 
principal value, with variables A and fj, restricted to the 
interval [—1,1], and the dependence on x and to is en- 
coded in the function 



V (A) = exp [iiruX + 7t VA 2 /2] 



(7) 



It is interesting to note that, after the substitution x 2 i— > 
—it, at <j) = 0, the integral kernel /C coincides with that 
arising in the calculation of time-dependent correlation 
functions of the one-dimensional interacting Bose gas at 
zero temperature ^5|. A comparison of the universal 
function c(x) as inferred from Eqs. 1(3} , (|5I7|I with the 
results of direct numerical simulation is shown in Fig. ^ 
Before outlining the derivation of these results, several 
remarks are in order: (i) The universality of Eqs. 10 can 
be inferred from the universality of Eq. (|13[) below [6|. 
As such, these results can be applied to the parametric 
evolution of spectra which obey Dyson statistics only "lo- 
cally" . (ii) Although we have not succeeded in obtaining 
a direct proof, in accordance with the conjecture made in 
Ref. [l3| , the distribution of single-level shifts appears to 
assume a Gaussian form p (lo, x) = / 2<J ( X ) / ^/2-na{x) 
at any value of x. The corresponding width of the Gaus- 
sian can be expressed as 

At small x, where p(u>,x) can be inferred from the level 

velocity distribution a ( x ) x2 j which reflects 

(identifying time with x 2 ) independent "diffusion" of in- 
dividual levels. In the opposite limit x t— > 00, making 
use of the known asymptotic dependence |(| c(x) 
— 1/tt 2 x 2 obtained from a perturbative analysis, one can 



The resulting strongly sub-diffusive behavior at large x 
can be ascribed to the rigidity of the spectrum "hemming 
in" the meandering levels, (iii) The generating function Z 
is a particular case of a more general object Z q (J, J'; x; cf>) 
which defines 

P q (J, J'; x) = P %^Z q (J, J'; x- <f>) (8) 

as the probability that the number n(J) of levels {ei(x)} 
in the (not necessarily contiguous) interval J and the 
number n'(J') of levels {ti(x + x)} in the interval J' dif- 
fer by exactly q. Taking J and J' to be semi-infinite 
intervals J e = (— 00, e] and J e + W = (— 00, e + w], one can 
derive generalizations of Eqs. Jjyl and (J3J which involve 
correlations between levels ti (x) and ei+ q (x + x) : 
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The exact analytic expression for Z will be given below, 
(iv) In some applications, the fixed perturbation XV may 
be of finite rank; i.e., it may possess only a finite num- 
ber r of non-zero eigenvalues In such situations 
Eqs. ((HI retain their validity providing T>(\) is replaced 
by e l7VujX det sc (l — iXTZ), where 1Z is the reactance ma- 
trix for scattering off the potential XV 13], and det sc 
denotes the determinant in the space of scattering chan- 
nels. In this case x may be identified with any of the 
variables parametrizing 1Z. (v) Despite the existence of 
well-developed analytical tools for the study of integral 
kernels with the structure of (|6a[) 15|, it is at present 
unclear whether these methods can be generalized to ac- 
commodate 4> integration in Eqs. and Q). 

The analysis of parametric statistics of individual en- 
ergy levels relies on a technical device which ensures that 
the level ei(x + x) is indeed the "descendant" of Ci(x) 
by demanding that it has the same ordinal number as 
counted from the bottom of the spectrum. Specifically, 
due to the absence of level crossings, the intra-level ve- 
locity correlation function coincides with the conditional 
average 
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c{x) = I du{5 n( j e) ^, {Je+w) 
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x d s 6 [e - e t (x)} d- x B [e + lo - e 3 {x + x)]), (11) 



where 8 n , n i denotes the Kronccker (5-symbol, and 9 is the 
step-function. The corresponding distribution of level 
shifts p{uo, x) is given by an analogous expression with d s 
replaced by d e (and no integration over lo). By generaliz- 
ing the corresponding non-parametric formula for P n (J) 
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[l8j| . our starting point is the general expression for the 
probability P nn r(J, J') to find n levels in the interval J 
of the unperturbed sequence and n! levels in the interval 
J' of the perturbed sequence, 
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Here r kk / represents the multi-point parametric correla- 
tion function of DoS 0, ^| integrated over the interval 
jk (g, j/k w j^ n ^ ne corresponding measures d/ij and dp, ji . 
Owing to the determinantal structure of the DoS corre- 
lation function, r kk i can be represented in the form of a 
fermionic functional integral 



r kk , = detK / X>4T>* 



J dvj(u)£(u)£(u)j (^J d(j,jt(w)ri(w)r)(w)J e 



f du f dWt>(u)K 1 (u,w)'S>(w) 
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where ^ = (£,??) is a fermionic doublet. Here 
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where the matrix elements of the operator sine kernel k 
of the unitary Dyson ensemble are k(u — w) = sin7r(u — 
w)/n(u - w), and [V Q k]{u,w) = l 2 ^ l du ^ k(u - w). 
Fixing the difference n' — n = q and summing over all n, 
one obtains the probability P q { J, J'; x) that the numbers 
of levels in the two intervals J and J' differ by q: 



P q (J,J';x) = J2 
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Substituting Eq. i fTSjl into Eq. one finds 



P q (J,J';x) = (-l) 
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and <Ji are the Pauli matrices. The determinants are 
understood as functional determinants on the space of 
two-component functions defined on the product interval 
J ® J' . The expression in curly brackets in Eq. I|ltj|) can 
be identified as the generating function Z . 

In order to apply Eq. I|ltj|) to the computation of the 
intra-level velocity correlation function and the distribu- 
tion of the parametric level shifts, one must set J = J e 



and J' — J t+UJ . Setting q — and using Eq. ijTTJl , one ob- 
tains Eqs. © and Q, where Z is identified with the first 
determinant in Eq. (|16f) . For q ^ 0, a similar procedure 
leads to Eqs. © and (|TU1) . 

The use of semi- infinite intervals to define n(J e ) and 
n'{J t +uj) is justified only if the support of the spectrum is 
finite. The latter condition would be trivially fulfilled if 
one were to use, instead of k(u — w), the exact Christoffel- 
Darboux kernel of the GUE whose scaling limits interpo- 
late between the sine kernel inside the Wigner semicircle, 
and the Airy kernel at its endpoints. However, in prac- 
tice, employing such a kernel would present significant 
technical difficulties. To circumvent this problem, we use 
a regularized kernel 

ks ( u - w )= siM ("- m ) e -dffl%l+H) i (17) 

7T (U — W) 

where the limit 5 — > is implied in all expressions in- 
volving this kernel. Using Eq. Ijl3(l it is easily shown that 
(n(Joo)) - (n'(Joo)) = 2/5, and ([n(Joo) - n'(Joc)] 2 } ~ 
0(6). Thus, although the regularization formally vio- 
lates the level number conservation, the corresponding 
error tends to zero in the limit 6 — > 0. In the following 
we will suppress the index 5. 

As written, Eq. i|16|) involves a matrix oscillating inte- 
gral kernel defined on a product of semi-infinite intervals. 
However, as we will now show for the case q = 0, it can 
be rewritten in the form of Eqs. © which is (i) more 
amenable to numerical analysis, and (ii) makes the in- 
tegrability of the kernel (in the sense discussed in Ref. 
|l5j ') manifest. Without loss of generality, we can set 
e = 0, and shift the variables so as to define the deter- 
minant on the quadrant (— oo,0] ® (— oo,0]. The corre- 
sponding shift operator is absorbed into the redefinition 
Vq — > V = e ujd ^ du T>o- The term involving the (^-function 
in the upper right corner of Eq. lfH|l can be separated to 
reveal the dyadic structure of the remainder: 
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Now, using the identities (Ico + zJD 1 lcr+) 1 = 

(1<t - Z-£>~ 1 ler+), and det(lcr + Z-X>~ 1 lo+) = 1, 
where a+ = (<Ji + i<T2)/2, one obtains 



Employing the Fourier representations, 



det 



lcr - KU (<j>) 



k{u) 
5{u) 



dX 



;i - ami + a) | eiXnu 



= det 1 - (z + fc - z+z-k \T> 1 t)v + z-(f) 1 k) V ) . and making use of Eq. 0, one finds 



z+zA (v^ij f> + z_ (i)' 1 ^ £>) (it, M)) 



-e 



i2 [ 27ri J-oo A — /i — iS 



d^—- -^V + z-V- 1 (X)V}. (18) 



Finally, the cyclic invariance of the determinant and the 
identity z + + z- = are used to perform the integrals 

in the u-w space, with the resulting kernel in the X-fi 
space having the form of Eqs. 10. Remarkably, taking 
the limit S — > in the final expressions leads to a non- 
singular kernel defined in terms of the Cauchy principal 
value integral Ij6b|l . 



As a final comment, it should be noted that the method 
of using the 4> integration to "count" the levels in con- 
junction with the regularization analogous to l|17|) is 
equally applicable to other Dyson ensembles. However, 
at present there exist no analogs of Eq. 1)1 -ill for other 
ensembles, and thus our consideration is perforce limited 
to the unitary case. 
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FIG. 1: Level velocity correlation function obtained from 
Eqs. 15I7H (solid line) vs. direct numerical simulation 
of large random matrices flrj| (dots). The width a (x) of 
the Gaussian distribution of level shifts together with the x 2 
asymptotics at small x is shown in the inset. 



